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ABSTRACT 

A  marriage  of  elastic-plastic  fracture  mechanics  techniques 
with  thermoplastic  finite  element  analyses  is  developed  to  examine 
crack  growth  in  the  presence  of  weld-induced  residual  stresses.  A 
hypothetical  crack  growth  relation  based  on  the  crack  tip  opening 
displacement  is  used.  Three  problem  areas  are  studied:  stress 
corrosion  cracking  in  a  girth-welded  pipe,  fatigue  crack  growth 
under  cyclic  loading  in  a  butt-welded  plate,  and  dynamic  crack 
propagation  under  Impact  loading  in  a  butt-welded  plate.  Compari¬ 
sons  with  computations  carried  out  under  conventional  linear  elas¬ 
tic  assumptions  are  made.  It  is  found  that,  in  all  three  cases, 
neglect  of  the  plastic  deformation  caused  by  the  welding  process 
appears  to  be  anti-conservative.  It  is  concluded  that  more  realis¬ 
tic  computations  for  crack  growth  in  and  around  welds  than  are  com¬ 
monly  used  may  be  needed  for  realistic  structural  Integrity 
assessments. 


INTRODUCTION 

A  significant  proportion  of  all  structural  failures  cm  be 
traced  to  cracks  emanating  in  and  around  welds.  Crack  growth  in 
welded  regions  must  be  strongly  affected  by  the  presence  of  the 
plastically  deformed  material  Indigenous  to  the  welding  process. 
Yet,  present  day  fracture  mechanics  analysis  procedures,  which  are 
largely  based  on  linear  elastic  conditions,  do  not  directly  treat 
such  complications.  While  elastic-plastic  fracture  mechanics  anal¬ 
ysis  procedures  have  been  developed,  they  have  previously  been 


applied  to  account  only  for  crack  tip  plasticity  itself.  Recently, 
a  further  step  has  been  taken  by  the  authors  through  the  use  of 
postulated  elastic-plastic  crack  growth  relations  for  crack  growth 
in  weld-induced  plastic  deformation  fields*-^,  This  work  is  sum¬ 
marised  and  assimilated  in  this  paper. 

The  immediate  objective  of  the  work  reported  in  References  1-3 
was  to  critically  examine  the  assumptions  of  linear  elastic  mate¬ 
rial  behavior  commonly  made  in  analysing  weld  cracking  problems. 
Three  separate  problems  were  addressed.  As  shown  in  Table  1,  these 
included  two  different  welded  structures  -  a  girth-welded  pipe  and 
a  butt-welded  plate  -  and  three  different  crack  growth  mechanisms  - 
stress  corrosion,  fatigue  and  unstable  crack  propagation.  The 
solution  procedure  and  the  individual  results  are  first  presented 
in  what  follows  with  general  conclusions  for  future  progress  in  the 
analysis  of  crack  growth  in  the  presence  of  weld-induced  residual 
stress  drawn  from  them. 


Table  1 .  Problems  Examined 


Structure 

Material 

Applied  Loading 

Simulated 

Cracking 

Mechanism 

Girth-welded 
pipe  with  circum¬ 
ferential  crack 
(axisymmetry) 

Type  304 

Stainless 

Steel 

Constant  Tension 

Stress  Corrosion 

Butt-welded 
with  edge  crack 
(Plane  Strain) 

HY-80  Steel 

Cyclic  Tension 

Fatigue 

Butt-welded  plate 
with  edge  crack 
(Plane  Strain) 

HY-80  Steel 

Impulsively 
Applied  Tension 

Dynamic  Fracture 

THE  ANALYSIS  PROCEDURE 

The  analysis  procedure  followed  in  this  work,  which  is  basi¬ 
cally  the  same  for  all  three  problems  examined,  consists  of  three 
main  steps: 


The  residual  stress  field  induced  in  a  welding  process  is 
computed  using  an  incremental  thermoplastic  finite  element 
analysis  procedure. 

2.  Crack  growth  is  simulated  by  sequential  node  release  along 
a  pre-set  crack  plana  located  in  the  weld  heat-affected 
cone. 

3.  A  postulated  elastic-plastic  crack  growth  relation  is  used 
to  infer  crack  length  as  a  function  of  time  or  loading 
history  to  simulate  a  particular  crack  mechanism. 

A  comparison  with  a  computational  result  made  using  a  commonly 
accepted  linear  elastic  approach  is  then  made  to  assess  the  signif¬ 
icance  of  the  linear  elastic  assumption  and  the  essential  neglect 
of  residual  plasticity  inherent  in  such  an  approach. 


Residual  Stress  Analysis 

The  residual  stress  analysis  procedure  is  one  that  has  been 
used  successfully  at  Battelle  in  a  variety  of  applications*-6. 

It  consists  of  two  parts.  First,  a  thermal  analysis  is  made  to 
obtain  the  time-temperature  history  for  each  point  in  the  body  for 
each  Individual  welding  pass.  Then,  these  histories  are  used  as 
input  to  an  incremental  elastic-plastic  finite  element  model  to 
determine  the  stress  and  deformation  state  of  the  weld  and  the  base 
material  as  the  weld  is  deposited.  Because  each  welding  pass  is 
considered  on  an  individual  basis,  the  residual  stress  and  strain 
state  that  exists  at  the  completion  of  one  weld  pass  constitutes 
the  initial  condition  for  the  next.  The  final  residual  stress 
state  is  that  which  exists  at  the  completion  of  all  of  the  weld 
passes. 

In  the  work  described  here,  the  finite  element  models  for  the 
welding  analysis  were  designed  with  a  line  of  double  nodes  along  a 
pre-set  crack  plane.  This  had  no  effect  on  the  residual  stress  dis¬ 
tribution  as  it  was  assumed  that  an  Initial  (small)  crack  appeared 
after  the  welding  process  was  complete'.  The  potential  crack  plane 
was  located  in  the  heat-affected  zone  ^ut  was  otherwise  arbitrary. 

The  two  welded  structures  examined  in  the  work  reported  here  - 
a  girth-welded  pipe  and  a  butt-welded  plate  -  are  shown  in  Figures 
1  and  2,  respectively.  The  corresponding  finite  element  models  are 
shown  in  Figures  3  and  4.  It  can  be  seen  that  a  number  of  simpli¬ 
fications  have  been  introduced  for  computation  convenience.  These 
Include: 


S  SECTION  OF  A  SEVEN-PASS  GIRTH-BUTT  WEID 
,  4-INCH  DIAMETER  SCHEDULE  80  TYPE  304 
NLESS  STEEL  PIPE 


CROSS  SECTION  OF  AN  18-PASS  BUTT  WELD 
IN  AN  HY-80  STEEL  SHIP  STRUCTURE 


•  weld  lint  symmetry 

t  tvo-dimtntiontl  dtforottion  (l.t.,  axisymaetry  or  pltnt 
•train) 

•  one-dimensional  crack  growth  (i.e.,  conctntrlc  or 
col linear). 

Note  that,  because  the  cracka  are  not  auppoaed  to  be  at  the  center 
of  the  weld  (cf,  Figurea  3  and  4),  the  flrat  of  theae  simplifying 
assumptions  means  that  two  cracks  exist  in  the  analysis. 

The  temperature-dependent  material  properties  used  in  the  anal** 
yses  of  the  pipe  weld  are  those  of  Type  304  stainless  steel  shown 
in  Figure  5.  The  properties  for  the  butt-welded  plate  are  those  of 
the  ship  structure  steel  HY-80  given  in  Figure  6.  In  the  first 
problem,  the  heat  Inputs  were  taken  from  an  actual  experience7 
while  in  the  second,  typical  values  were  used.  These  together  with 
the  weld  pass  and  structure  geometry,  suffice  to  determine  the 
residual  stress  distribution.  The  normal  atresses  acting  on  the 
prospective  crack  plane  for  the  two  welded  configurations  are  shown 
in  Figures  7  and  8.  It  can  be  seen  that,  while  these  are  quite 
different,  they  share  a  high  tensile  stress  near  the  surface. 

Thus ,  an  edge  crack  would  be  likely  to  grow,  particularly  when  the 
residual  stresses  are  abetted  by  a  tensile  applied  stress. 


Subcritical  Crack  Growth  Analysis 

The  equivalence  between  the  crack  tip  opening  displacement 
(CTOD)  and  the  stress  Intensity  factor  in  small-acale  yielding  is 
well  known.  Recent  progress  in  elastic-plastic  fracture  mechanics 
has  further  revealed  the  distinctive  role  played  by  the  crack  tip 
crack  opening  displacement  in  crack  initiation  and  stable  growth 
in  large  scale  yielding  conditions®"^-®.  Specifically,  for  the 
initiation  of  crack  growth,  the  CTOD  can  be  expressed  as 


r 

EY 

small-scale  yielding 

u 

deformation  plasticity 

where  K  is  the  stress  intensity  factor,  J  is  the  J-integral  param¬ 
eter,  E  is  the  elastic  modulus,  Y  is  the  yie.’.d  stress  while  a  and 
dn  are  numerical  constants  on  the  order  of  unity.  In  addition, 
for  extended  stable  crack  growth,  the  CTOD  appears  to  take  on  a 
constant  value.  While  certainly  not  conclusive  evidence  that  the 
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FIGURE  8.  RESIDUAL  STRESS  CN  POTENTIAL  CRACK  LINE  IN  A  1-INCH  THICK 
BUTT-WELDED  PLATE  OF  HY-80  STEEL 


Nominal  Stress 


CTOD  is  the  controlling  parameter  for  subcritical  crack  growth  as 
well,  for  lack  of  an  alternative,  it  will  be  so  taken  in  this  work. 
Note  that,  because  of  the  equivalence  represented  by  Equation  (1), 
this  choice  is  not  Inferior  to  one  based  on  either  K  or  J  in  any 
event . 


Starting  from  an  assumed  initial  crack,  crack  growth  is  simu¬ 
lated  in  each  of  the  finite  element  models  by  sequential  node  re¬ 
lease  along  the  line  of  double-noded  elements.  Each  nodr.  pair  is 
released  by  diminishing  the  Initial  force  that  exists  between  them 
to  zero.  This  is  done  over  from  five  to  ten  load  increments.  The 
value  of  6  for  a  given  crack  length  is  then  the  value  of  the  CTOD 
that  exists  when  the  load  has  vanished.  Hence,  it  is  a  value  at 
one  finite  element  spacing  behind  the  actual  crack  tip.  The  re¬ 
sults  obtained  from  the  finite  element  models  are  shown  in  Figures 
9  and  10.  Note  that  in  the  latter  problem  two  solutions  are  shown: 
one  with  no  applied  stress,  the  second  with  an  applied  tensile 
stress  normal  to  the  crack  plane  of  67  percent  of  the  room  tempera¬ 
ture  yield  stress  of  HY-80  steel  (80  ksi). 

Figures  9  and  10  contrast  the  CTOD  values  obtained  by  advanc¬ 
ing  the  crack  through  the  finite  element  model  under  two  different 
conditions.  First,  the  computed  weld-induced  plastic  deformation 
is  left  unaltered  and  an  incremental  plasticity  computation  made. 
This  is  the  elastic-plastic  analysis.  Second,  a  simplified  ap¬ 
proach  i6  followed  wherein  (1)  only  the  normal  component  of  the 
residual  stress  acting  on  the  potential  crack  plane  is  retained, 
and  (2)  linear  elastic  behavior  is  assumed.  This  is  denoted  as  the 
simple  elastic  analysis  and  typifies  that  commonly  used  for  this 
kind  of  problem. 

Stress  corrosion  cracking  is  often  supposed  to  occur  according 
to  a  power  law  relation  of  the  type 


da 

dt 


CK® 


(2) 


where  a  denotes  the  crack  length,  t  is  time,  while  C  and  m  are 
material  constants.  Because  this  relation  is  clearly  valid  only 
under  small  scale  yielding  conditions,  the  relation  between  K  and  6 
expressed  by  Equation  (1)  can  be  used  to  cast  it  into  the  equiva¬ 
lent  form 


da 

Ji 


C’6b/2 


(3) 


where  C'  is  also  a  material  constant*.  Having  6*5  (a)  from 
Figure  9,  Equation  (3)  can  be  integrated  numerically  to  find  the 
crack  length  in  terms  of  a  dimensionless  time  t* 


CTOD,  8,  ft  in. 


COMPUTED  CRACK  TIP  OPENING  DISPLACEMENT  IN  A  BUTT-WELDED  HY-80  PLATE  AS  A  FUNCTION  OF  CRACK  GROWTH 
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•  .m/2-1 
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o 

where  h  denotes  the  wall  thickness.  The  result  is  shown  in 
Figure  1 1 . 

Notice  that  the  results  given  in  Figure  11  were  obtained  with¬ 
out  the  imposition  of  an  applied  stress.  Thus,  despite  the  fact 
that  the  residual  stresses  are  self-equilibrating,  catastrophic 
crack  growth  can  occur.  The  reason  is  that  the  stresses  are  redis¬ 
tributed  as  crack  growth  occurs  and,  at  least  until  a  net  compres¬ 
sive  force  is  freed  by  the  growing  crack,  a  positive  crack  driving 
force  exists;  cf,  Figure  9.  Of  more  Importance,  the  results  of 
Figure  11  indicate  that  the  simple  elastic  analysis  predicts  an 
anti-conservative  result  in  that  the  time-to-failure  If  greater 
than  in  the  more  rigorous  elastic-plastic  analysis. 

Fatigue  crack  growth  under  a  uniform  cyclic  loading  can  often 
be  adequately  characterized  in  the  form 

§  ■ c  <K».*  -  W  <» 

where  N  denotes  a  load  cycle  number  and  C  and  m  are  material 
constants.  Again  introducing  the  CTOD  from  equation  (1)  gives 

where  6xu*x  *n^  5 min  *re  the  CTOD  values  that  would  be  attained 
under  the  maximum  and  minimum  load  levels,  respectively.  Conse¬ 
quently,  the  number  of  cycles  required  to  achieve  a  given  crack 
length  can  be  found  by  integrating  Equation  (6)  via 


•o  <61/2  -  61/2) 
v  max  min' 

The  results,  using  the  CTOD  values  given  in  Figure  10,  are  shown  in 
Figure  12.  It  can  be  seen  that  the  simple  elastic  analysis  is  once 
again  anti-conservative. 


Elastic  -  Plastic 
Analysis 


BUTT-WELDED  PLATE 


The  governing  relatione  for  unstable  crack  propagation  and 
arrest  in  elastodynamic  conditions  are 

K  -  Kn  ,  a  >  0 

D  (8) 

K  <  KD  ,  a  -  0 


where  Kq  la  known  as  the  dynamic  fracture  toughness.  Once  again, 
so  long  as  small  scale  yielding  conditions  are  satiafled,  such 
relations  are  equivalent  to  thoae  couched  in  terms  of  the  CTOD.  In 
particular,  the  critical  CTOD  value  for  dynamic  crack  propagation 
would  then  be  simply 


6 


D 


(9) 


Clearly,  if  a  time  varying  loading  is  imposed  on  a  cracked 
body,  the  CTOD  value  will  also  vary.  So  long  as  6  <  {q,  the 
crack  tip  will  be  stationary.  If  6  becomes  equal  to  £q,  crack 
growth  will  occur  at  a  rate  such  that  the  equality  is  maintained. 
Arrest  occurs  when  the  equality  can  no  longer  be  satisfied. 


A  dynamic  computation  was  performed  in  which  the  ship  struc¬ 
ture  shown  in  Figure  2  was  subjected  to  a  suddenly  imposed  load  of 
42  ksi,  about  SO  percent  of  the  room  temperature  yield  stress. 

This  load  was  held  for  30  Msec  and  then  dropped  to  aero.  A  value 
of  Kc  •  37.9  ksiVin .  was  used  to  reflect  the  lower  toughness 
existing  in  the  heat-affected  cone.  Using  Equation  (9),  this  gave 
a  critical  CTOD  value  of  0.0006  inch.  Both  a  linear  elasto¬ 
dynamlc  and  an  elastic-plastic  dynamic  calculation  were  made.  The 
elastic-plastic  analysis  was  based  on  the  entire  residual  stress 
and  deformation  field  together  with  Incremental  dynamic  plasticity. 
However,  in  accord  with  common  practice,  the  residual  stresses  were 
completely  ignored  in  the  elastic  analysis.  The  computed  CTOD 
values  for  the  two  analyses  are  shown  in  Figure  13. 


It  can  be  seen  in  Figure  13  that  the  elastic-plastic  analysis 
predicts  that  unstable  crack  growth  would  quickly  occur  (at  approx¬ 
imately  15  Msec).  It  also  predicts  that  the  crack  would  penetrate 
the  wall.  In  contrast,  the  elastic  analysis  does  not  predict 
initiation  of  growth  until  much  later  (about  40  Msec)  and  predicts 
arrest  soon  thereafter.  Consequently,  the  simpler  procedure  is 
once  again  found  to  be  anti-conservative. 


DISCUSSION 


It  is  Important  to  recognize  that  the  prime  purpoae  of  this 
work  waa  not  to  arrive  at  quantitative  reaults  for  different  types 
of  material  behavior.  Rather,  it  waa  to  critically  examine  a  set 
of  assumptions  that  are  commonly  used  in  analyzing  a  class  of  prob¬ 
lems  for  given  material  behavior.  The  specific  materials  consid¬ 
ered  in  this  study  are  very  ductile  and  tough  (e.g.,  Type  304 
stainless  steel,  HY-80  steel)  which  undoubtedly  exacerbates  the 
differences  that  ware  found.  It  is  quite  possible  that  other 
materials— and,  perhaps  more  importantly,  residual  stress  fields 
Induced  without  large-scale  plaatic  deformations— would  show  con¬ 
siderably  less  difference.  Of  equal  Importance,  since  the  elastic- 
plastic  crack  growth  criteria  needed  for  the  purposes  of  this  study 
have  not  been  established,  a  pragmatic  approach  was  taken  to  obtain 
comparative  results.  All  of  these  factors  would  be  borne  in  mind 
in  interpreting  the  results  given  in  this  paper. 

The  basic  assumption  that  has  been  called  into  question  here 
is  the  applicability  of  linear  elastic  fracture  mechanics  in  the 
presence  of  weld-induced  residual  stress  fields.  This  has  been 
addressed  by  performing  two  parallel  computations  where  this 
assumption  has  and  has  not  been  made.  Hence,  while  there  are 
undeniably  many  aspects  of  the  calculations  that  can  be  improved 
upon,  because  the  two  computations  were  otherwise  performed  on 
exactly  the  same  basis,  these  cannot  be  of  critical  importance. 
Indeed,  the  comparison  has  revealed  such  wide  disparities,  that, 
neglect  of  the  Inelastic  deformation  accompanying  welding  would 
appear  to  be  unequivocally  incorrect. 

Conversely,  the  work  presented  in  this  paper  should  not  be 
taken  as  a  blanket  indictment  of  LEFM-based  crack  growth  predic¬ 
tions.  The  mathematical  convenience  of  LEFM  is  too  useful  to  not 
play  an  Important  part  in  the  assessment  of  weld  cracking.  What 
appears  to  be  needed  is  some  sort  of  plasticity-enhanced  LEFM 
procedure,  possibly  calibrated  with  the  more  rigorous  analyses 
described  in  this  paper,  that  can  be  confidently  applied  even  in 
the  presence  of  large-scale  plasticity  and  its  attendant  residual 
stress  fields.  How  this  can  best  be  done  is  an  open  question  at 
this  time.  However,  as  the  work  reported  herein  so  strongly  sug¬ 
gests,  the  necessity  for  it  is  not. 


CONCLUSIONS 

Elastic-plastic  fracture  mechanics  research  has  identified  the 
CTOD  as  a  key  crack  growth  parameter.  The  use  of  this  finding  in 
conjunction  with  thermoplastic  finite  element  analysis  procedures 
has  enabled  more  realistic  computations  of  crack  growth  in  the 


presence  of  weld-induced  residual  deformation  and  stresses  to  be 
made.  Comparisons  of  these  results  with  commonly  used  approaches 
based  on  linear  elastic  fracture  mechanics  indicate  that  the  latter 
could  be  highly  anti-conservative. 
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